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Abstract 



> 

This paper is a study of some of the critical properties of a simple model for 
' flux. The model is motivated by gauge theory and is equivalent to the Ising 

Q\ ' model in three dimensions. The phase with condensed flux is studied. This 



is the ordered phase of the Ising model and the high temperature, deconfined 



phase of the gauge theory. The flux picture will be used in this phase. Near 



the transition, the density is low enough so that flux variables remain use- 



^ ' ful. There is a finite density of finite flux clusters on both sides of the phase 



transition. In the deconfined phase, there is also an infinite, percolating net- 
work of flux with a density that vanishes as T — > T+. On both sides of the 
critical point, the nonanalyticity in the total flux density is characterized by 
the exponent (1 — a). The main result of this paper is a calculation of the 
critical exponent for the percolating network. The exponent for the density 
of the percolating cluster is £ = (1 — a) — ((p — 1). The specific heat exponent 
a and the crossover exponent cp can be computed in the e-expansion. Since 
C < (1 — a), the variation in the separate densities is much more rapid than 
that of the total. Flux is moving from the infinite cluster to the finite clusters 
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much more rapidly than the total density is decreasing. 
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I. INTRODUCTION 



This paper is a study of some of the critical properties of a simple model for flux. The 
model is motivated by gauge theory and is equivalent to the Ising model in three dimensions. 
The main emphasis is on the phase in which the flux is condensed. This is the ordered phase 
of the Ising model and the high temperature, deconfmed phase of the gauge theory. In such 
a phase, there is a percolating network of flux in addition to clusters of finite size. 

The character expansion for the Ising model is a product over links. Each term in the 
product has the form (1 + zss'). The expansion of the product can be viewed as a sum over 
configurations of links that have the zss' term rather than the 1. By definition, these are 
the links with flux. This machinery is commonly used to calculate expansions in the small z 
phase. I will use it to describe configurations in the large z phase also. Naturally, it cannot 
be the basis for a perturbative expansion in that phase. In the large z phase, there is a 
percolating network of links with flux. This network must disappear as z — > z+. The main 
result of this paper is a calculation of the critical exponent for the vanishing density of the 
percolating network. 

All of the calculations are for the model introduced above. My interest in the model 
originates in the subject of finite-temperature gauge theory. In the low-temperature, confined 
phase of a non-Abelian gauge theory, it is believed that sources are connected by a tube of 
color-electric flux rather than by a dipole Coulomb field. This flux-tube description of the 
low-temperature phase has been used fLj to account for the properties of the second-order, 
finite-temperature, deconfining phase transition of pure-glue, SU(2), gauge theory. At high 
temperature, a plasma of free gluons with perturbative corrections can account for some 
properties. However, there are also nonperturbative effects in the deconfmed phase. The 
correct physical picture for T — > T c + is not established. I will use the flux picture near 
the phase transition in the deconfmed phase. In this phase, the flux is condensed, but the 
density may be low enough so that flux variables remain useful. I expect a finite density of 
flux clusters of a finite size on both sides of the phase transition. In the deconfmed phase, 
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there is also a percolating network of flux with a density that vanishes as T — > T c + . 

In the strong-coupling description of SU(2) lattice gauge theory, the flux on links is 
labeled by the irreducible representations of SU(2). At each site, the representations must 
combine to have a singlet piece. For SU(2), this can happen only if there is an even number 
of links with flux at each site. 

These configurations of flux on links are subject to both quantum and finite-temperature 
fluctuations. The quantum fluctuations reflect the fact that these configurations are eigen- 
states of the Hamiltonian only in the lowest order of the strong coupling expansion. The 
magnetic term in the Hamiltonian moves the flux on the lattice. Thus, the flux is thickened 
into tubes. In the usual, but unproven, picture of confinement, that is as far as it goes. The 
flux is not spread into Coulomb-like configurations, and tubes of flux remain a useful de- 
scription. On a scale larger than the basic length 1/Aqcd, the eigenstates of the gauge field 
Hamiltonian are configurations of flux tubes, and the energy is proportional to the length 
of the flux. There is an additional assumption that on these large scales, it is sufficient to 
include only the basic flux strength associated with fundamental sources. 

Finite-temperature fluctuations are the main interest of this work. As the temperature 
increases, small loops and other clusters of flux appear with low density. Increasing temper- 
ature brings increasing density. Most work fl[ with this kind of picture has concentrated on 
the region below the critical temperature. At the critical temperature, there is a condensa- 
tion of flux. I assume that the flux picture can be carried into the high-temperature phase 
where the condensation appears as a percolating network of flux. This paper deals with the 
region just into the high-temperature phase. 

With this background, the simplified model that will be studied can be presented. Space 
is a three-dimensional cubic lattice. There are variables 9 on links that can have the values 
1 or to indicate the presence or absence of flux. The energy of a link with flux is a. The 
weight for a configuration is 
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The sum over configurations is restricted to those in which the number of links at a site that 
have flux is even. Let the collection of all such configurations be C It is a subset of the 
unrestricted collection of all configurations C . The partition function is 

Z = Je-fE^W. (1.2) 

c 

Some ideas from percolation theory will be used to study this model. They are introduced 
in the next section. Also, it will be shown that the flux model is equivalent to the Ising 
model. In the third section, the main result is the exponent for the vanishing density of the 
percolating flux as T — > T+. 

Poo ~{T-T c f (1.3) 

with 

C=(l-a)-fo>-l). (1.4) 

The specific heat exponent a and the crossover exponent ip can be computed in the 
e-expansion. The fourth section contains a brief summary. 

II. PERCOLATING FLUX 
A. General remarks 

The discussion of the previous section has led to a model for flux defined by the partition 
function 

Z = Y,e~^ ae{l) . (2.1) 

a 

The link variable 6(1) = 1 or to indicate the presence or absence of flux on link I. The 
sum is restricted to configurations in which the number of links at a site that have flux is 
even. In the high-temperature phase, there is a condensation of flux. If the density of the 
condensed flux is not too high, then a description in terms of flux variables will be useful. 
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To discuss this phase, I will draw on ideas from percolation theory 0. With- 
out the restriction on configurations, the model becomes pure bond percolation with 
p = e~ a l T j{\ + e~ a l T ) as the probability for flux on a link. In percolation, both phases 
have a finite density of flux clusters of finite size. In the large p phase, there is also a single 
percolating cluster. Since p is the probability that a link has flux, p is also the total density 
of flux per link. Thus, in terms of the control parameter p, the total density has exponent 1, 
i-e. ptotai = Pc + (p — Pc) 1 - The critical value of p c is about 0.25, so the critical link density is 
also about 0.25. Since just a quarter of the links have flux, this remains a useful description 
of configurations. 

The density of the infinite cluster vanishes with a rate given by the exponent f3 ~ 0.403. 
It is important to note that the exponent (3 is smaller than one. This means that there must 
be compensating behavior in the density of finite clusters for p > p c . My interest is in the 
analogous relations for the simplified flux model. 

The flux model differs from percolation only through the restriction on configurations. 
To study the flux model, the ideas from pure percolation can be applied. In particular, 
assume that there is a single percolating cluster in the high temperature phase and that the 
nonanalyticities at the phase transition are characterized by critical exponents. Naturally, 
the numerical values of the exponents will be different from those of percolation. 

B. Equivalence with the Ising model 

Now it is time to show that the model for flux is equivalent to the Ising model. To do 
this, introduce site variables to enforce the restriction on configurations, and then do the 9 
sums. The site variables will turn out to be the Ising spins. 

First, consider the sum of 9(1) over the Z's contained in the set of links L(i) that hit a 
site i: 

E m (2.2) 

l€L{i) 
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To force this to be even at each site, introduce the site variables s(i), which take the values 
±1. Then the factor 

I £ s(z)^'^ m (2.3) 

s(i)=±l 

has the desired effect. A factor like this is introduced into the partition function sum for 
each site i. Then the restriction on the 9 sums can be relaxed. With the abbreviations 

£ = II(!£)> (2-4) 

s i s(i) 

z = e~ a l T , and ss' for the product of the pair of spins on the sites that bound a given link, 
the partition function is 



z = e e n e ~ Tam n s«^ 6L(!) 9{l) (2.5) 

SCI i 



i 

= EEIK*«O' (0 (2.6) 

s C I 

=En( i +^ / ). (2.7) 

s I 

This is a version of the Ising model 

Y^e^ 88 ' (2.8) 

s 

with z = tanh (3. The links with flux are those that receive the term zss' rather than 1 in 
the expansion of Y\i(l + zss'). 

C. Additional remarks on SU{2) gauge theory and the Ising model 

Thus, the simple model for SU(2) flux is equivalent to the Ising model. The appearance 
of the Ising model is not surprising because the gauge theory and the Ising model are related 
by another route 0. Finite-temperature, SU(2), gauge theory has a global Z(2) symmetry. 
The order parameter for the second-order, gauge theory phase transition is the Wilson 
line. The effective theory for the lines inherits the Z(2) symmetry. This effective theory 
is in the universality class of the three-dimensional Ising model. An approximate, strong 



coupling, character expansion for the effective theory involves a similar, but somewhat more 
complicated, appearance of flux on links. 

In Ref. [|j], there are discussions of the adjoint Wilson line that involve the behavior of 
flux. The finite-temperature expectation value of the adjoint line has the form || 

(Li) = c + d ± \^\ 1 - a . (2.9) 

In the deconfined phase, this receives contributions from the finite clusters and from the 
percolating cluster of flux. However, the separate contributions need not have the exponent 
(I — a). This paper will address that issue. To the extent that the simple model captures 
the essential aspects of the gauge theory, then the conclusion will be that the two gauge 
exponents are also different. 

Now that the general picture has been outlined, the task of the next section is to deter- 
mine the behavior of the flux around the phase transition. The main emphasis will be on 
the critical exponent for the density of the percolating cluster. 



III. CALCULATIONS 

The purpose of this section is to calculate the critical exponents for the total flux density 
and for the density of the percolating cluster. The model has been formulated in terms of 
flux on links 

Z = Y j e-^ am (3.1) 

c 

and as the Ising model 

^zna+^o- (3.2) 

The discussion begins with the second form. Operators associated with the desired flux 
densities must be identified. In order to separate the percolating cluster, the spins will be 
generalized to unit vectors s that can point along either direction of iV axes. Then the 
iV — > 1 limit will be taken to return to the original case. This is a technical device that is 
closely related to a similar trick in ||. 

8 



A. Total flux density 



The total flux density can be obtained straightforwardly without the generalization of 
variables mentioned above. Consider a particular link, the spin variable pair ss f on the sites 
that bound it, and the operator (z — ss'). In the product over all links, this will multiply the 
factor (1 + zss') for the link in question to give (z 2 — l)ss'. This eliminates configurations 
without flux on this link. Configurations with flux on this link receive a factor of (z 2 — 1) 
rather than z in their weight. Therefore, the expectation value of the operator 

X = —^—(z-ss>) (3.3) 
[z 2 - 1) 

gives the probability that a given link has flux. 

Evidently, X is closely related to the internal energy per link E = —ss'. Near the critical 
point 

- E = c + d±\z- z c \ x - a (3.4) 

with c pa 0.332, z c pa 0.218, and a pa 0.12 0. Thus, the operator X that gives the total flux 
density has the critical value X c = z c (z c — c)/(z 2 — 1) = 0.0261 and the critical exponent 
(l-a) w 0.88. 

These numbers can be compared with analogs from pure percolation. At the critical 
point, the factor for an additional link of flux is z c pa 0.22. This is slightly larger than 
the analogous number 0.21 for a self-avoiding walk || and is somewhat smaller than the 
factor p c / (1 — p c ) pa 0.33 in pure percolation. The critical density of flux in the Ising 
model is about one tenth of what it is in pure percolation. Since many clusters that are 
possible in percolation are not possible in the Ising model, the total Ising density should be 
lower. For example the smallest cluster in percolation has just one link, while in the Ising 
model, it has four. With the same factor for an occupied link, this will give a much more 
dilute concentration of occupied links in the Ising model. In the growth of a large cluster, 
the excluded volume will be smaller in the Ising model so that the entropy per link of a 
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large cluster will be larger in the Ising model. Since criticality occurs when the entropy 
and energy balance, the critical value of the link factor will be smaller in the Ising model 
than in percolation. This further reduces the critical density in the Ising model relative to 
percolation. 

This concludes the initial discussion of the total flux density. 

B. Percolating flux density 

In the high temperature phase, the percolating cluster and the finite clusters contribute 
to the total density. The goal is to determine the exponent for the percolating cluster near 
the critical point where its density goes to zero. It will turn out to be a bit smaller than the 
(1 — a) for the total density. It follows that there must be a complementary increase in the 
contribution from the finite clusters. Thus, both components vary more rapidly than their 
sum near the critical point. More flux is moving from the percolating cluster to the finite 
clusters than is disappearing. 

The density of the percolating cluster can be separated by using a somewhat indirect 
method. Generalize the spins to unit vectors s that are restricted to point in one of the 2N 
directions that are parallel or antiparallel to the axes of an iV-dimensional space. Generalize 
the partition function to 

^Ena+^-s')- (3.5) 

s I 

Now Y, s means sum each spin over the 2N allowed orientations. This defines the generalized 
spin model. 

As before, the links with flux are those that have the term with the zss' rather than 1 in 
the expansion of the product. All the spins in a cluster must be along the same axis. Thus, 
there are N mutually repelling flavors of Ising flux. 

In the z > z c , condensed phase, an infinitesimal external field will determine a preferred 
flavor and break the cubic symmetry. I assume that there is a single percolating cluster on 
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which the spins are restricted to the preferred flavor. To be specific, let the preferred flavor 
be the first spin component. 

To measure the density of the percolating cluster, consider the spin s at a given site. It 
satisfies the simple identity 

S 2 = E^ = 1- (3-6) 

a 

From the first spin component si, construct the operator 

Y = s{ - i. (3.7) 

Its average will be zero in the symmetric phase. In the percolating phase, the finite clusters 
are still symmetric, so they do not contribute to (Y). On sites hit by the percolating cluster, 
Y = (N — 1)/N. Thus, 

(Y) = x (site density of the percolating cluster) (3.8) 

There are a number of equivalent forms for Y: 

Y = si - i (3.9) 

= 4-^4 (3-io) 

-^si-^^sl (3.11) 

(3.12) 



N 1 N Z^°a 



The last equality is true if we limit ourselves to matrix elements of Y with other operators 
that contain si only. This is sufficient for our purposes. The expectation value of the 
operator Y = (sf — s\) gives the density of the percolating cluster. 
Recall that the goal is to find the exponent C, in 

(F)~(T-T C ) ? . (3.13) 

For (T — T c ) < where there is no percolating cluster, (Y) = 0. The exponent is related to 
the anomalous dimension of the operator in the usual manner 
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C = udy (3.14) 

= u(d-2 + j Y ) (3.15) 
= (l-a)-(p-l). (3.16) 



The combination = v(2 — jy) is a crossover exponent |10|. It will turn out that cp > 1 



Thus, the decrease in the density of the percolating cluster is much more rapid than that of 
the total density. 

Note that the operator Y gives the density of sites hit by the percolating cluster. Earlier, 
it was link densities that appeared. The ratio of sites hit to links hit is certainly between 
1/3 and 1. Also, the density of the percolating cluster approaches zero. It follows that the 
critical exponents for the percolating site and link densities must be the same. For the finite 
clusters, where the density does not approach zero, the site and link exponents might be 
different. 

The task at hand is to calculate the anomalous dimension of Y and then take N 1. 
To do this, I will pass to a field theory description and use the e-expansion. 

C. Field theory 

The next step is to find a field theory version of the generalized spin model in which the 
anomalous dimension of Y can be computed. First, ideas from the renormalization group 
and universality will be used as the basis for a guess. Then a systematic transformation and 
approximation will be done to show that the guess can be put on a stronger foundation. 

A field theory in the universality class of the spin model should be on a three-dimensional, 
Euclidean space and should have a field with the internal symmetry of the order parameter. 
The generalized model has iV-dimensional cubic symmetry. A simple guess is iV-component, 
scalar field theory with "cubic anisotropy" |IT]. The Lagrangian density is 



C = Ud<P) 2 + >V + W) 2 + & E# (3-17) 
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We also need an operator to play the role of Y. An obvious combination to try is <p\ 



This can be related to the more general and much-studied [K| 

M N 



B = ^[(N-M)J2<f>l-M £ <&] 

■ 1 a=M+l 



(3.18) 



a=l 



We are interested in the specialization M — > 1 followed by A" — > 1 at the end. If we stick to 
matrix elements with <f>\, then 



B 



(3.19) 



Besides the evident relation to Y and Y, the operator in (|3.18| ) has satisfying technical 
properties. It does not mix with </> 2 under renormalization, so that it has its own anomalous 
dimension. Also, its anomalous dimension is independent of M. 

To put this guess in a more respectable position, a generalized Gaussian transformation 



12] can be done. First, the partition function for the generalized model must be put in the 



form of a sum over configurations with a Boltzmann weight e~ u l T . Then a conventional 
path to a field theory can be followed. Since the model does not satisfy (s ■ s') 2 = 1, the 
Hamiltonian must be more complicated than H cx J2i s ■ s' . The two-term form 



H/T = J2[as- s' + b(s- s 



is sufficient. For the partition function, this gives 



A2l 



(3.20) 



s I 



-[as-s'+b{s-s') 2 ] 



(3.21) 



If this is to agree with ( |3.5| ), we need 



z = — tanh a and e = cosh a. 



(3.22) 



For the Gaussian transformation, a slightly different expression for summand in ( |3.20| ) 
is a better starting point flT3"||: 



(3.23) 
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Now follow the method described in Ref . |12j . This involves the introduction of two pairs of 
iV-component, real- valued fields. The pair u a and 4> a is associated with s a , and the pair v a 
and Xa is associated with s 2 a . In the end, there will be a field theory for alone. 
First Z is rewritten as 

Z = jd(j) Jd X Jdu JdvY.e-^^'^^ J]eS a { ^ W[Ua(i) - Sa(i)]+XaWK(i) - s ' (i)]} , (3.24) 

s % 

then the sum on the spin variables s is done. This gives 

fd X Jdu /d we -E ! (^'+^)+Ej^»-(i)+x(0-«] 

(IIEXI! ^ U) ) [II 2 ^sh <p ak (k)]. (3. 25) 

i o-i j k 

The next step is to do the x an d v integrals. These are in a factor that is a function of the 
indices {aj} 

Z x = Jd X y^e-Ei^'+EjxW-^-x^Wl (3.26) 
= e -& (3.27) 

The is the Kronecker delta. The subscripts i and i' refer to the sites that bound 

the link /. Notice that this is the Boltzmann weight for an iV-state Potts model with site 
variables aj. This could also have been obtained directly from ( |3.23| ). Next include the cij 
sums and the cosh factors to get a functional of (j>. 

Z * = (nE) e ^ E;5(ai ' ai,) [Il2cosh0 afc (fc)] (3.28) 

i a-i k 

= M(II 2cosh <MA0» (3-29) 

k 

The quantities Z P and ((•)) refer to the partition function and expectation value for the Potts 
model. 

Note that this is an antiferromagnetic Potts model. It is in the disordered phase. Equa- 
tion ( |3.22| ) shows that, when z and a are small, b ~ a 2 /2. So b will be much smaller than 
a, and the second term in ( |3.2(J| ) is a small correction that can be computed in a strong- 



coupling series. Without the second term in (|3.20| ), the mean field estimate is a c ~ N/6. 
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This gives b « 0.13 for JV = 3, and b ta 0.056 for N = 2. For A" = 1, the actual a c « 0.22 
can be used to give b ~ 0.024. These are small numbers. For comparison, consider that the 
critical values for the antiferromagnetic Potts model are b c ~ 0.82 for N = 3 and b c ~ 0.44 
for N = 2 [|n|]. The N = 1 case is trickier. The N — 1, ferromagnetic Potts model is 
percolation. The critical probability 1 — e bc = p c « 0.25 corresponds to b c ~ —0.29. The 
N = 1, antiferromagnetic case is percolation with negative probability. The power series in 
p that gives p c ~ 0.25 has all positive coefficients through order p 10 . However, an analysis 
of the series indicates the presence of a weak singularity for negative p that is a little closer 
to the origin ]nj . There is no evidence for a singularity that would correspond to b as small 
as 0.024. 

The contribution that Z 2 makes to the action for will be computed through order 4 . 
In the end, it will be easy to see that higher powers of are irrelevant. 

After expanding to order 4 , computing the Potts expectation values, and exponentiat- 
ing, the 0-dependent part of the exponent is 

-u = + &<(<)]» - H¥IM) 2 } + IEE«3^(0J<(7)»e- (3.30) 

The subscript c on the last Potts average indicates that it is connected. Using the fact that 
all values of are equally likely in the symmetric phase, the first three terms become 

HK(») + - §[£ E 5#(01 2 }- (3-31) 

i a a 

The summand of the last term in —U becomes 

iE^O'Mi) (3-32) 

where 

= «<5(a, ai)«J(6, 6j)» c (3.33) 

is the connected Potts two-point function. 

With the Potts model away from its critical point, D is short-ranged. Since it is the 
long-wavelength part of that determines its critical behavior, we can take to be constant 
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over the range of \i — j\ for which D is significantly different from zero. The first correction 
is an irrelevant operator with spatial derivatives and four powers of (p. That leaves 

EiE^W^WIE^j) ( 3 - 34 ) 

i ab jj^i 

for the last term in —U. The internal index structure of D a b has the symmetry of the 
symmetric phase of the Potts model. This allows for a term that is independent of the 
indices a and b and for a term that is proportional to S(a, b). Thus, the terms that appear 
in (|3.34f) are of the same form as those in ( |3.31|) . The only effect of the term b(s ■ s') 2 is to 



alter the coefficients of the 4 terms in the potential. 

Now we see that the contribution from U to the action for cf) is a 4 potential with cubic 
anisotropy. Returning to ( |3.25|) , we have 

Z = (dcj) Jdue~^i aM+ ^ ${l)<i) - u . (3.35) 

The last step is to do the u integral. However, there is little point in carrying the calculation 
further. It is clear that, except for the 6-dependence of the coefficients in U, this effective 
theory for is the same as that which would have come from the conventional, bilinear 
form a(s ■ §'). It is a textbook result []E| that the u integral and another long-wavelength 
approximation give a "kinetic" term for and a <fi 2 term to add to U. 

The conclusion of this discussion is that a Gaussian transformation confirms the original 
guess. The field theory associated with the generalized spin model is an iV-component, 4 , 



scalar field theory with cubic anisotropy. Since this is a well-studied example [11], we can 
draw from established results. 



D. Fixed points and exponents 

Now that the preliminaries are out of the way, we can get to the essentials relatively 
easily. The heavy lifting has already been done by others and is described rather nicely in 
]nj . For this field theory, there are five possibilities revealed by the e-expansion. There can 
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be a fluctuation-induced, first-order phase transition. Since the Ising and the SU (2) phase 
transitions are both second order, this case can be eliminated. The second possibility is the 
Gaussian fixed point. This has free field exponents and is unstable in both directions in the 
coupling constant plane. Both these properties make it unlikely to be the case of interest. 
The bare action of Q3.17 ) begins away from the origin of the coupling plane. The Ising and 



gauge models have nontrivial fixed points. 

There are three nontrivial fixed points called Ising, cubic, and Heisenberg. At the Ising 
fixed point, the N components decouple into N independent fields each at the 1-component 

fixed point. There is one stable and one unstable direction. Again, it is unlikely that 
( 3.17p will hit it. It has already been noted that the Ising flavors of the generalized spin 



model repel each other. Equivalently, ( |3.17|) includes coupling that is unstable at this fixed 



point. However, ( |3.17|) is an approximation, and there could be some fine cancellation that 
has been spoiled. If this were the correct fixed point, then the anomalous dimensions for 
<f) 2 and B would be the same, and the exponents for the percolating and the finite clusters 
would be the same. 

The cubic fixed point has the original cubic symmetry, while the Heisenberg fixed point 
has full O(N) symmetry. The relative stability of these two depends upon the value of N. 
For our case with N — > 1, the Heisenberg point is stable, while the cubic one has one stable 
and one unstable direction. Again, in the absence of some special cancellations that the 
approximations have missed, the cubic fixed point is not the one of interest. Also, in the 
iV — > 1 limit, this fixed point joins the Gaussian fixed point. The anomalous dimension for 
(f) 2 is zero. Since this is not the correct behavior for the total density, this cannot be the 
correct fixed point. 

Finally, we come to the stable and, therefore, most likely case: the Heisenberg fixed 
point. The e-expansions for the anomalous dimensions of <fi 2 and B are in Ref. ||10|| . The 
lowest order results for N — > 1 lead to 

a = e/6 and (p-l = e/18. (3.36) 
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Since ( — (1 — a) — (</? — 1), this gives the exponent for the percolating flux. For e = 1, the 
result is C = (1 — a) — 1/18. This may not be particularly precise. However, it reveals the 
important qualitative fact that ( is less than (1 — a). 

IV. SUMMARY 

In the low temperature phase, there are only finite clusters. In the high temperature 
phase, the total flux is the sum of the finite clusters and the percolating cluster. On both 
sides of the critical point, the nonanalyticity in the total flux density is characterized by the 
exponent (1 — a): 

Ptotai = c + d ± \T-T c \ 1 - a . (4.1) 

The main result of this work is the exponent for the vanishing density of the percolating 
flux as T -> T+. 

Poo ~ (T - T c f (4.2) 

with 

C=(l-a)-(<p-l). (4.3) 

The specific heat exponent a and the crossover exponent ip can be computed in the 
e-expansion. 

Since ( < (1 — a), the variation in the separate densities is much more rapid than that of 
the total. Flux is moving from the infinite cluster to the finite clusters much more rapidly 
than the total density is decreasing. 
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